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ABSTRACT 

Potential  genetic  differences  between  healthy  and  declining  sugar  maple 
{Acer  saccharum  Marsh.)  were  investigated  at  four  sites  in  Ontario.  Eleven 
polymorphic  enzyme  loci  were  assayed  using  starch-gel  electrophoresis  and  the 
genetic  data  was  analyzed  statistically  in  conjunction  with  decline  inde.x  (DI)  ratings 
of  individual  trees.  Significant  allelic  heterogeneity  between  healthy  and  declining 
subsamples  was  observed  at  one  locus  (Mdh-3)  at  one  site  (Davies),  but  there  was 
no  evidence  of  direct  involvement  of  Mdh-3  in  resistance  or  susceptibility  to 
decline.  Levels  of  expected  heterozygosity  of  healthy  and  declining  subsamples  were 
also  significantly  different  at  the  Davies  site  with  the  declining  subsample 
demonstrating  greater  heterozygosity.  Results  from  the  other  sites  were 
contradictory  but  not  statistically  significant.  No  trend  was  apparent  concerning 
observed  heterozygosity.  Comparisons  of  allele  frequencies  and  genetic  diversities 
among  healthy  and  declining  subsamples  of  an  individual  site  may  be  confounded 
by  cohort  associations  or  family  structure.  Therefore,  it  is  desirable  to  limit  such 
structure  when  sampling,  but  this  may  be  difficult  in  natural  populations  where 
parentage  is  unknown.  Fixation  indices  (f  js)  were  consistently  greater  in  declining 
subsamples  than  in  the  healthy  subsample  of  the  same  site  suggesting  that  the 
proportion  of  heterozygous  loci  per  individual  relative  to  Hardy-Weinberg 
expectations  tends  to  be  greater  in  the  healthy  individuals.  Such  a  relationship 
could  occur  if  inbreeding  leads  to  a  predisposition  to  decline  and/or  if  individual 
heterozygosity  contributes  to  resistance  to  decline.  The  relative  contribution  of 
either  of  these  possible  conditions  to  the  observed  trend  in  F i^  is  not  easily 
discerned. 
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INTRODUCTION 

Over  the  past  few  years  forest  decline  has  attracted  global  concern.  Largely 
enigmatic  in  origin,  it  has  been  described  throughout  much  of  North  America  and 
Europe  (Hornbeck  and  Smith,  1985;  Knight,  1987;  Linzon,  1985;  Schiitt  and 
Cowling,  1985).  The  changing  chemical  climate  has  received  considerable  attention 
as  a  probable  contributing  factor.  It  is  generally  agreed  that  air  pollutants  such  as 
ozone  and  oxides  of  nitrogen  and  sulphur  accompanied  by  an  increase  in  acidity 
of  precipitation  are  detrimental  to  the  world's  forests  (Anderson,  1989;  Linzon, 
1985;  Morrison,  1984;  Schulze,  1989).  Although  acidic  deposition  has  often  been 
singled  out  as  a  primary  malefactor,  it  is  unlikely  that  it  alone  is  the  major  cause 
of  forest  decline  (Pitelka  and  Raynal,  1989;  Roberts  et  al.,  1989;  Skelly,  1989). 
Decline  may  be  the  result  of  the  interactions  of  a  number  of  stresses  of  which 
pollution  may  be  one  part  (Bardy,  1989;  Linzon,  1986;  Mueller-Dombois,  1988; 
Pierre  and  Queiroz,  1988).  In  some  species,  decline  may  be  a  stage  in  a  normal 
process  initiated  by  periodic  natural  disturbances  such  as  severe  drought  (Hosking 
and  Hutcheson,  1988)  or  may  be  a  response  to  long  term  climatic  shifts  (Hamburg 
and  Cogbill,  1988).  Whatever  its  cause,  the  problem  of  forest  decline  is  real  and 
widespread. 

Among  the  Canadian  species  evoking  the  greatest  concern  is  sugar  maple 
(Acer  saccharum  Marsh.).  Reports  of  maple  decline  in  Ontario  began  as  early  as 
1947  and  many  of  those  which  have  followed  have  listed  its  cause  as  unknown 
(Mcllveen  et  al.,  1986).  It  is  not  surprising  that  maple  syrup  producers  have 
blamed  acid  rain  for  the  decline  of  their  sugar  bushes  since  acidic  deposition  in 
southern  Ontario  is  among  the  greatest  in  North  America  (Wisniewski  and  Kinsman, 
1988),  but  indications  are  that  acidic  precipitation  is  but  one  more  stress  in 
addition  to  those  of  insect  outbreaks  and  drought  (Hendershot  and  Jones.  1989; 
Linzon,  1986;  McLaughlin  et  al.,  1985).  According  to  Hendershot  and  Jones  (1989) 
areas  with  poor  soils  and  marginal  climate  are  predisposed  to  decline. 
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Previous  research  in  forest  trees,  which  has  primarily  considered  air 
pollution  to  be  the  causal  agent,  has  pro%ided  reason  to  believe  that  genetic 
differences  may  exist  between  healthy  and  deteriorating  individuals.  Variation  in 
tolerance  levels  to  pollutants  has  been  recorded  among  provenances  in  Pinus 
strobus,  (Kriebel  and  Leben,  1981),  P.  sylvestris  (Huttunen.  1978).  and  Picea  abies 
(Huttunen,  1978;  Scholz  and  Venne.  1989),  and  among  clones  in  Picea  abies  (Scholz 
and  Venne,  1989),  Pinus  sircbus  (Houston  and  Stairs,  1973),  and  Populus 
tremuloides  (Wang  et  al.,  1986).  Concern  about  the  potential  erosion  of  genetic 
resources  through  the  elimination  of  sensitive  individuals  has  been  expressed 
(Karnosky  et  al.,  1989;  Kreibel  and  Leben.  1981;  Scholz.  1981).  The  purpose  of 
the  present  study  was  to  determine  whether  possible  genetic  differences  between 
healthy  and  declining  sugar  maple  within  the  same  sites  are  discernable  at  the 
allozyme  level.  Regardless  of  the  specific  causal  agents,  the  determination  of  any 
link  between  genetic  markers  and  visible  injury  could  provide  a  foundation  for 
selection  and  management  decisions. 

MATERIALS  AND  METHODS 

Daca  Collection 

The  sugar  maple  examined  in  this  study  were  from  four  sites  (Figure  1) 
that  had  been  considered  in  an  earlier  study  (Beak  Consultants  Ltd.,  1989).  .At 
each  site,  bud  tissue  was  collected  from  24  sugar  maple  trees  previously  classed  as 
healthy  and  from  24  trees  (23  at  Boothby)  which  had  been  classed  as  declining, 
giving  a  total  sample  size  of  47  to  48  trees  per  site.  Bud  tissue  was  prepared,  and 
starch-gel  electrophoresis  was  carried  out  as  previously  described  (Perry  and 
Knowles,  1989).  A  total  of  nine  enzyme  systems  were  assayed  (Table  1).  Only 
variable  loci  (eleven  in  total)  were  considered  and  although  Sod-2  was  found  to 
have  low  level  variability  in  earlier  studies  (Leiliot,    1987;   Perry  and  Knowles, 


1989),  it  was  invariant  in  the  presently  sampled  trees.  Therefore  it  was  excluded 
from  further  analysis. 


FIGURE  1.  LOCATIONS  OF  FOUR  SUGAR  MAPLE  SITES  SAMPLED  IN  ONTARIO. 


TABLE  1^  ENZYME  SYSTEMS  ASSAYED  IN  BUD  TISSUE  OF  SUGAR  MAPLE. 


ENZYME  SYSTEM 


Aspartate  aminotransferase  (AAT) 
Aldolase  (ALD) 
Diaphorase  (DIA) 
Leucine-amino  peptidase  (LAP) 
Malate  dehydrogenase  (MDH) 
6-Phosphogluconic  dehydrogenase  (6PG) 
Phosphoglucose  isomerase  (PGI) 
Phosphoglucomutase  (PGM) 
Superoxide  dismutase  (SOD) 


E.G. 

BUFFER 

DESIGNATION 

SYSTEM^ 

2.6.1.1 

B 

4.1.2.13 

H 

1.6.4.3 

B 

3.4.11.1 

B 

1.1.1.37 

H 

1.1.1.44 

H 

5.3.1.9 

B 

2.7.5.1 

H 

1.15.1.1 

B 

Buffer  systems  as  described  by  Cheliak  and  Pitel  (1984). 

A  decline  index  (DI)  was  available  for  each  sampled  tree.  This  index  was 
based  on  the  percentage  of  dead  branches  in  the  tree  crown,  the  percentage  of 
slight  and/or  strong  chlorosis,  and  the  percentage  of  undersized  leaves  (Mcllveen 
and  McLaughlin,  1988).  The  DI  values  adopted  were  the  averages  of  estimates 
made  by  two  observers. 

Statistical  Analysis 

Allele  frequencies  were  estimated  for  both  subsamples  of  each  population 
based  on  the  observed  genotypic  arrays.  In  many  cases,  small  expected  numbers 
(<5)  precluded  valid  testing  of  allelic  heterogeneity  using  standard  chi-square  or 
C-tests.  Therefore,  heterogeneity  among  allele  frequency  estimates  within  sites  was 
evaluated  via  Monte  Carlo  methods  as  follows.  Using  a  computer  program,  for  each 
site  and  each  locus,  a  genotypic  distribution  equivalent  to  the  pooled  observed 
healthy  and  declining  samples  was  randomly  divided  into  two  synthetic  classes.  The 
proportion  of  10,000  such  synthetic  divisions  having  allele  frequencies  equally  or 
more  unevenly  distributed  than  those  observed  was  considered  as  the  probability 
of  the  observed  distribution. 


Expected  heterozygosities  (//g)  ^^r^  estimated  for  each  locus  within  each 
subsample  as: 

where  p^  was  the  estimated  frequency  of  the  /th  allele  at  the  locus.  A  measure  of 
the  deviation  of  observed  heterozygote  frequency  from  that  expected  under 
Hardy-Weinberg  equilibrium,  F[s  was  calculated  for  each  subsample  as: 

where  //q  was  the  observed  proportion  of  heterozygotes.  The  contribution  of  each 
locus  to  a  subsample's  F[s  was  weighted  by  its  expected  heterozygosity  in  that 
subsample.  Values  of  f[s  greater  than  0  indicate  an  excess  of  homozygous 
genotypes  relative  to  Hardy-Weinberg  expectations  while  negative  values  are 
indicative  of  an  excess  of  heterozygotes.  For  each  site,  subsample  differences  in 
//e,  Hq,  and  f  IS  were  evaluated  using  paired  /-tests. 

Some  overlap  of  DI  between  subsamples  within  three  of  the  four  sites  was 
noted.  For  example,  the  healthy  subsample  at  the  Boothby  site  had  DI  values 
ranging  from  1.87  to  17.39  while  the  declining  subsample  had  trees  with  DI 
estimates  ranging  from  12.14  to  44.84.  To  remove  possible  effects  of  this  overlap, 
all  trees  within  a  site  were  ordered  according  to  their  DI  and  the  twenty  trees  with 
the  highest  DI  within  a  site  were  considered  as  a  new  declining  subsample  and  the 
twenty  trees  with  the  lowest  DI  within  that  site  were  reclassified  as  a  new  healthy 
subsample.  Herein  these  new  subsamples  are  referred  to  as  reduced  data.  The 
statistical  analyses  described  above  were  performed  both  on  these  reduced  data 
and  on  the  data  of  subsamples  as  originally  designated  (original  data). 

The  possibility  of  a  relationship  between  heterozygosity  and  vigour  was 
explored  through  a  regression  analysis  of  individual  tree  data.  The  independent 
variable  was  the  proportion  of  loci  heterozygous  in  a  tree  and  the  dependent 
variable  was  the  tree's  DI.  This  analysis  was  performed  for  each  site  using  the 
computer  program  MCREG  of  BIOITAT  I  (Pimentel  and  Smith,   1985). 


RESULTS 

Observed  allele  frequencies  of  the  healthy  and  declining  subsamples  (original 
data)    within    the    four    sites    are    presented    in  Table  2.    Monte  Carlo  methods 


TABLE  2.  ALLELE  FREQUENCIES  OF  HEALTHY  (H)  AND  DECLINING  (D)  SUGAR 
MAPLE  AT  FOUR  SITES  (ORIGINAL  DATA). 


ALLELE 

site' 

BOTH 

BR.A.K 

DAV 

LMF 

LOCUS 

H 

D 

H 

D 

H 

D 

H 

D 

Aat-1 

0.98 

0.94 

0.96 

0.96 

1.00 

0.92 

1.00 

0.96 

0.02 

0.02 

0.04 

0.04 

0.00 

0.06 

0.00 

0.02 

0.00 

0.04 

0.00 

0.00 

0.00 

0.00 

o.oo 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.00 

0.00 

Aat-2 

0.69 

0.85 

0.63 

0.73 

0.90 

0.79 

0.79 

0.83 

0.31 

0.15 

0.31 

0.25 

0.08 

0.21 

0.21 

0.17 

0.00 

0.00 

0.06 

0.02 

0.02 

0.00 

0.00 

0.00 

Ald-3 

0.90 

0.89 

0  98 

0.90 

1.00 

0.92 

0.90 

0.94 

0.10 

0.09 

0.02 

0.10 

0.00 

0.08 

0.08 

0.06 

0.00 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.00 

Dia-S 

0.8S 

0.96 

0.96 

0.92 

0.96 

0.92 

0.85 

0.92 

0.08 

0.04 

0.04 

0.08 

0.04 

0.08 

0.13 

0.06 

0.04 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.02 

Lap-1 

0.98 

0.98 

1.00 

0.94 

1.00 

1.00 

1.00 

1.00 

0.02 

0.02 

0.00 

0.06 

0  00 

000 

0.00 

0.00 

Lap-2 

0.96 

1.00 

0.98 

0.98 

1.00 

1.00 

1.00 

0.98 

0.04 

0.00 

0.02 

0.02 

0.00 

0.00 

0.00 

0.02 

Mdh-3 

0.92 

1.00 

0.94 

0.96 

0.98 

0.85* 

0.96 

0.96 

0.08 

0.00 

006 

0.04 

002 

0  15 

0.04 

002 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

6pg-l 

1.00 

1.00 

0.98 

1.00 

1.00 

0.98 

0.96 

1.00 

0.00 

0.00 

0.02 

0.00 

0.00 

0.02 

0.04 

0.00 

Pgi-2 

0.88 

0.76 

0.73 

0.79 

0.73 

0.67 

0.85 

0.79 

0.06 

0.17 

0.21 

0.19 

008 

0.15 

0.13 

0.10 

0.04 

0.04 

0.02 

0.00 

0.15 

0.06 

0.02 

004 

0.00 

0.00 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.02 

0.02 

0.02 

0.04 

0.12 

0.00 

0.06 

Pgm-1 

0.90 

0.91 

0.96 

0.90 

0.81 

0.88 

0.90 

0.88 

0.00 

0.02 

0.02 

0.02 

0.02 

0.04 

0.04 

0.00 

0.10 

0.07 

0.02 

0.08 

0.17 

0.08 

0.06 

0.12 

Pgiii-2 

0.98 

1.00 

1.00 

0.98 

1.00 

0.98 

1.00 

0.96 

2 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.00 

0.02 

3 

0.02 

0.00 

0  00 

0  02 

000 

000 

0  00 

002 

^Site  names:  BOTH=Boothby;  BRAK  =  Brackenbridge,  DAV=Davies;  LMF  =  L.M.  Frost. 
Significant  (P<0.05)  difference  between  subsample  allele  frequencies  within  site. 
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demonstrated  that  heterogeneity  was  significant  (P=0.04)  at  Mdh-3  in  the  Davies 
site.  Allele  frequencies  in  the  reduced  data  (not  presented)  varied  little  from  those 
of  the  original  data  and  conclusions  regarding  heterogeneity  were  unaltered. 

Mean  DI  of  subsamples  within  sites  ranged  from  5.72  to  30.98  in  the 
original  data.  When  the  reduced  data  were  viewed,  this  range  increased  from  4.74 
to  33.55.  In  general,  declining  subsamples  had  mean  DI  about  4  times  greater  than 
the  corresponding  healthy  subsamples  using  original  data  (Table  3).  When  reduced 
data  were  considered,  the  DI  of  declining  subsamples  were  just  over  5  times  greater 
on  average. 


TABLE  3.  MEAN  DAMAGE  INDEX  OF 
HEALTHY  AND  DECLINING  SUGAR  MAPLE 
AT  FOUR  SITES  (ORIGINAL  DATA). 


SITE  HEALTHY        DECLINING 


BOOTHBY  7.40  25.54 

BRACKENBRIDGE  5.72  27.45 

DAVIES  8.55  30.98 

L.M.  FROST  6.72  29.86 

POOLED  SITES  7.10  28.49 


Paired  /-tests  indicated  a  significant  difference  in  the  levels  of  e.xpected 
heterozygosity  of  the  healthy  and  declining  subsamples  within  the  Davies  site.  The 
declining  subsample  had  a  greater  H^  in  both  the  original  and  reduced  data  (Table 
4).  However,  in  the  remaining  sites  H^  was  usually  larger,  but  not  significantly 
so,  in  the  healthy  subsamples.  No  consistent  trend  was  noted  concerning  a  possible 
relationship  between  DI  and  observed  heterozygosity  and  none  of  the  sites 
demonstrated  significant  heterogeneity  of  Hq  among  subsamples  (Table  5). 
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TABLE  4.  EXPECTED  HETEROZYGOSITIES  PER  LOCUS  (Hg)  OF  HEALTHY  AND 
DECLINING  SUGAR  MAPLE  AT  FOUR  SITES. 


ORIGINAL  DATA 

REDUCED  DATA 

SITE 

HEALTHY 

DECLINING 

HEALTHY 

DECLINING 

BOOTHBY 

0.17 

0.13 

0.18 

0.14 

BRACKENB  RIDGE 

0.16 

0.17 

0.17 

0.16 

DAVIES 

0.12 

0.19 

0.12 

0.19 

L.M.  FROST 

0.16 

0.15 

0.15 

0.15 

POOLED  DATA 

0  16 

0  16 

0  16 

0.16 

Estimates  differ  significantly  (P<0.05)  between  subsamples  within  site. 


TABLE  5.  OBSERVED  HETEROZYGOSITIES  PER  LOCUS  (Hq)  OF  HEALTHY  AND 
DECLINING  SUGAR  .MAPLE  AT  FOUR  SITES. 


ORIGINAL  DATA 

REDUCED  DATA 

SITE 

HEALTHY 

DECLINING 

HEALTHY 

DECLINING 

BOOTHBY 
BRACKENBRIDGE 
DAVIES 
L.M.  FROST 

0.18 
0.18 
0.11 
0.17 

0.14 
0.16 
0.17 
0.16 

0.19 
0.19 
0.12 
0.15 

0.14 
0.15 
0.17 
0.15 

POOLED  DATA 


Regardless  of  the  direction  of  the  shift  in  H^  or  H ^  among  subsamples 
within  a  site,  f  jg  was  consistently  greater  in  the  declining  subsamples  than  in  the 
corresponding  healthy  subsamples  (Table  6).  In  addition,  when  compared  to 
estimates  obtained  using  original  data,  in  three  of  the  four  sites  Fjs  of  healthy 
subsamples  was  smaller  (more  negative)  and  that  of  declining  subsamples  was  larger 
when  made  with  reduced  data,  indicating  a  possible  strengthening  of  the  trend 
when  the  overlap  within  the  original  data  was  removed.  Although  paired  r-tests  did 
not  find  any  of  the  observed  differences  in  f  is  within  sites  to  be  statistically 
significant,  the  apparent  trend  was  compelling. 
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TABLE  6.  FIXATION  INDICES  (Fjg)  OF  HEALTHY  AND  DECLINING  SUGAR  MAPLE  AT 
FOUR  SITES. 


ORIGINAL  DATA 

REDUCED  DATA 

SITE 

HEALTHY 

DECLINING 

HEALTHY 

DECLINING 

BOOTHBY 

-0.054 

-0.036 

-0.083 

-0.020 

BRACKENBRIDGE 

-0.130 

0.075 

-0.169 

0.098 

DAVIES 

0.041 

0.095 

-0.008 

0.101 

L.M.  FROST 

-0.066 

-0.012 

-0.028 

-0.020 

POOLED  DATA 

-0.022 

0,050 

-0.044 

0.063 

Not  unexpectedly,  the  results  of  the  regression  analysis  reflected  the 
observations  made  of  H^  above.  The  possibility  of  a  relationship  between 
heterozygosity  and  state  of  decline  was  strongest  in  the  Davies  site  with  a 
correlation  (r)  between  individual  heterozygosity  and  DI  of  0.27  (P=0.06).  The 
Boothby,  Brackenbridge  and  L.M.  Frost  sites  had  /•  values  of  -0.22  (P=0.13),  -0.21 
(P=0.i6),  and  0.05  (P=0.74)  respectively.  Clearly  there  was  no  consistent  trend. 


DISCUSSION 

Several  observations  made  in  the  course  of  this  study  are  of  interest.  First, 
where  others  report  increased  heterozygosity  in  forest  trees  tolerant  of 
environmental  stress  (Geburek  el  al.  1987;  Miiller-Starck,  1985,  1989)  we  found  no 
evidence  of  such  a  trend  either  in  comparison  of  mean  heterozygosity  measures  of 
subsampies  or  in  the  regression  analyses.  In  contrast,  the  only  statistically 
significant  relationship  observed  was  contradictory  with  a  greater  H^  found  in  a 
the  declining  subsampie  of  the  Davies  site.  But  this  trend  was  not  observed  in  data 
of  the  other  sites.  It  should  be  noted  that  not  all  researchers  have  found  a  clear 
relationship  between  heterozygosity  and  apparent  tolerance  to  en\  ironmental  stress 
in  forest  trees  (Bergmann  and  Scholz,   1989). 
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The  Davies  site  was  also  unique  in  that  significant  allelic  heterogeneity  was 
observed  among  its  subsamples  while  no  such  heterogeneity  was  detected  for  the 
other  sites.  In  other  forest  trees,  heterogeneity  of  allele  frequencies  has  been 
observed  among  groups  sensitive  or  tolerant  to  environmental  stress  (Bergmann 
and  Scholz,  1985,  1989;  Geburek  et  al.,  1987;  MiJller-Starck,  1985,  1989;  Scholz  and 
Bergmann,  1984)  and  the  researchers  involved  have  often  suggested  that  the 
specific  enzymes  may  be  involved  in  resistance  mechanisms  or  are  acted  upon 
through  genie  selection.  We  have  no  evidence  that  the  Mdh-3  locus  is  directly 
involved  in  resistance  or  susceptibility  to  sugar  maple  decline.  Note  that  the 
distribution  of  allele  2  at  Mdh-3  in  the  Boothby  site  is  opposite  to  that  observed 
at  the  Davies  site.  Therefore,  if  the  decline  occurring  at  both  sites  is  attributable 
to  the  same  stress  factors,  Mdh-3  cannot  be  directly  involved  with  a  resistance 
mechanism  or  in  some  other  manner  be  a  specific  target  of  selection. 

Differences  in  allele  frequencies  between  tolerant  and  sensitive  phenotypes 
in  this  and  similar  allozyme  investigations  could  be  accounted  for  by  influences 
other  than  selection  imposed  by  environmental  stress.  We  suspect  that  the  allelic 
heterogeneity  observed  between  healthy  and  declining  subsamples  at  the  Da\  ies  site 
could  be  a  reflection  of  the  population's  genetic  architecture.  Sugar  maple  stands 
may  be  composed  of  a  number  of  cohorts  (progeny  produced  during  the  same  year) 
which  may  be  genetically  differentiated  (Mulcahy,  1975).  Willis  and  CotTman 
(1975)  suggest  that  stands  may  be  spatially  contagious  with  respect  to  age.  In 
addition,  through  allozyme  analysis  of  three  northern  sugar  maple  stands,  we  have 
found  evidence  of  spatial  genetic  structure  that  we  interpret  to  be  a  reflection  oi 
cohorts  or  family  groups  (manuscript  in  preparation).  Therefore,  if  cohorts  or 
family  groups  that  differ  in  their  allozyme  characteristics  also  differ  in  their 
resistance  to  pollution  or  other  stress  factors,  significant  differences  in  allozyme 
frequencies  or  diversities  between  tolerant  and  sensitive  groups  may  result.  These 
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allozyme  differences  would  likely  not  be  the  result  of  direct  action  of 
environmental  stress,  but  would  be  artifacts  of  cohort  association. 

Allelic  heterogeneity  observed  among  tolerant  and  sensitive  subsets  in  some 
earlier  studies  might  be  explained  in  a  similar  manner.  For  example,  where  seed 
is  collected  from  only  a  small  number  of  trees  and  the  resulting  progeny  are  tested 
for  resistance  to  a  particular  agent  (eg.,  SO2),  it  would  not  be  unexpected  that  the 
tolerant  subset  would  be  dominated  by  members  of  a  few  specific  families  and  the 
sensitive  subset  by  members  of  a  few  others,  if  resistance  is  a  heritable  trait.  At 
the  same  time,  allozyme  differences  are  expected  among  families  although  these 
differences  need  not  be  associated  with  resistance.  Clearly  it  is  important  to  limit 
family  structure  within  test  populations  in  order  to  avoid  possible  confounding 
effects.  This  is  difficult  in  natural  populations  where  parentage  is  unknown. 

Some  researchers  have  expressed  concern  that  some  alleles  are  only 
represented  in  the  sensitive  subpopulations  that  they  have  studied  and  thus  are 
considered  endangered  (Bergmann  and  Scholz,  1985,  1987,  1989;  Geburek  et  al., 
1987).  Considering  the  sample  sizes,  the  present  study  does  not  provide  evidence 
that  any  of  the  observed  alleles  are  at  risk.  While  several  alleles  were  found 
exclusively  in  the  declining  subsamples,  their  frequencies  were  low.  Actual  overall 
occurrence  of  the  most  frequent  of  such  subsample  specific  alleles  (allele  3  at 
Aat-1)  was  in  three  heterozygotes  out  of  the  total  sample  of  191.  This  is  not 
convincing  evidence  of  an  allelic  association  with  sensitivity  and  imminent  loss  of 
an  allele  since  the  probability  that  all  three  of  the  observed  heterozygotes  would 
be  grouped  together  in  the  same  decline  index  class  by  chance  alone  is 
approximately  0.12. 

What  inferences  can  be  drawn  from  our  findings?  Perhaps  the  most  notable 
trend  observed  was  that  of  F,s.  That  ail  four  populations  displayed  greater  f  ,5  in 
declining  subsamples  (Table  6)  may  not  be  convincing  on  its  own.  However,  when 
subsample  composition  was  redefined  in  accordance  with  D!  (reduced  data)  there 
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was  a  concurrent  decrease  of  f  ;s  in  healthy  subsamples  and  increase  of  F^^  in 
decliaing  subsamples  for  three  of  the  four  populations.  We  regard  this  as  evidence 
that  the  perceived  trend  is  something  more  than  an  ephemeral  aberration. 

This  trend  in  Fjs  indicates  that  the  proportion  of  heterozygous  loci  per 
individual  relative  to  Hardy-Weinberg  expectations  tends  to  be  greater  in  the 
healthy  individuals.  Such  a  condition  could  occur  if  the  more  inbred  individuals 
are  predisposed  to  decline,  or  if  heterozygosity  in  general  confers  resistance  to 
decline,  or  both.  It  is  difficult  to  distinguish  between  these  possibilities  (Mitten  and 
Jeffers,  1989).  That  the  significant  difference  in  heterozygosities  observed  among 
subsamples  at  the  Davies  site  had  the  declining  subsample  with  the  highest 
heterozygosity  would  suggest  that  heterozygosity  per  se  does  not  necessarily  impart 
resistance.  However,  if  all  inbred  individuals  were  eliminated,  the  lowest  fig  value 
expected  would  be  0  unless  other  factors  were  involved.  Lower  F^^  values  in 
healthy  or  tolerant  subsamples  might  tend  to  lead  to  higher  heterozygosities  in  those 
subsamples  than  in  declining  or  sensitive  subsamples,  but  not  always  as  evidenced 
by  the  Davies  site.  If  one  were  to  interpolate  these  observations  concerning  f  ,5  to 
a  breeding  program  with  a  goal  of  reducing  susceptibility  to  decline,  one  might 
select  superior  parents  based  on  field  observations,  and  then  make  selections  within 
families  according  to  allozyme  heterozygosity.  However,  it  should  be  noted  that 
further  study  is  required  to  confirm  the  proposed  trend. 

A  question  that  remained  following  an  earlier  study  (Perry  and  Knowles. 
1989)  was  whether  the  levels  of  variation  observed  in  five  sugar  maple  populations 
at  the  northern  range  limit  in  Ontario  was  representative  of  the  species.  Although 
the  data  of  the  present  study  were  not  collected  specifically  to  address  this  question 
and  the  sampling  methods  used  could  introduce  some  bias,  we  can  make  some 
tentative  observations.  In  the  northern  populations  previously  studied.  H^  ranged 
from  0.098  to  0.132  per  population  with  an  overall  mean  of  0.110.  In  the  present 
study,  if  the  same  loci  which  include  five  monomorphic  loci  (Mdh-1  and  2.  Sod-1, 


16 

2,  and  3)  are  considered,  H^  ranges  from  0.103  to  0.128  with  an  overall  mean  of 
0.115  (data  not  presented).  Therefore,  there  is  no  evidence  that  the  marginal 
populations  are  less  variable  than  central  ones  as  has  been  suggested  by  studies  in 
other  forest  tree  species  (Yeh  and  Layton,   1979;  Guries  and  Ledig,  1982). 

This  investigation  raises  many  new  questions.  Is  there  variation  in  the  level 
of  decline  among  cohorts  within  sites?  Obviously  tree  ages  would  have  to  be 
considered  to  address  this  question.  Are  declining  and  healthy  individuals 
distributed  randomly  within  sites?  This  would  require  that  tree  positions  be  mapped 
so  that  statistical  procedures  such  as  spatial  autocorrelation  analysis  could  be 
applied.  If  the  decline  has  a  spatial  element,  is  it  related  to  genetic  or 
environmental  factors?  How  might  artifacts  of  family  structure  or  cohort 
association  be  minimized  when  comparing  declining  and  healthy  subsamples  from 
natural  stands?  As  is  the  case  with  the  study  of  interactions  of  environmental  stress 
and  genetic  components  in  general,  the  problems  involved  in  the  investigation  of 
the  association  of  genetic  parameters  and  the  maple  dieback  are  multifarious. 
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